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Abstract 

Robust open-loop steering of a finite-dimensional quantum system 
is a central problem in a growing number of applications of information 
engineering. In the present paper, we reformulate the problem in the 
classical control-theoretic setting, and provide a precise definition of 
robustness of the control strategy. We then discuss and compare some 
significant problems from NMR in the light of the given definition. 
We obtain quantitative results that are consistent with the qualitative 
ones available in the physics literature. 

1 Introduction 

We consider an isolated n-dimensional quantum system with time evo- 
lution described by the following Schrodinger equation: 

ih\m) = Ham. (i) 
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Here \ip(t)) is a vector of unit norm in C n representing the state of 
the system at time t. The unitary time evolution of the system is 
governed by the system Hamiltonian: 

m 

H(t) = H + Y,H jUj (9,t). (2) 

3=1 

The internal Hamiltonian Hq £ C nx " is an Hermitian matrix describ- 
ing the free evolution of the system. The control Hamiltonian 

m 

H c (t) = Y,H jUj (9,t), 
i=i 

where Hj £ c nxn are a i so Hermitian matrices, accounts for the ef- 
fects of the control inputs ui(0,t), ...,u m (6,t) on the dynamics of the 
system. We assume that these control functions depend on a finite 
number of real parameters 9 = (9%, ...,8 p ), 9k £ T, with T being an 
open set in W. This kind of assumption is reasonable if we think of 
the small set of parameters we can control in an experimental setting. 

We consider the problem of steering the system from a given initial 
state |V>o) = 1^(0)) to a final state where |V>o) and are unit 
vectors in C n . We assume that the transition occurs (at t=T) when 
9 = 6* which we take as "nominal" value of the parameters. Clearly, 
if 9 9*, the transition will, in general, not occur. 

It is convenient to introduce the error probability for each control 
strategy. Consider the normalized final state for the time evolution, 
\tp(T, 9)}. It can be written as \ip(T,9)) = (tp 1 \ip(T,9))\ip 1 ) + \^ ± (9,T)) 
with \ifi (9, T)} orthogonal to If we imagine to perform a discrete 
measure 1 on an observable that has as eigenstate, the probability 
to obtain the eigenvalue associated |^i)(that corresponds to the prob- 
ability of finding the system in immediately after the measure) 
is: Pu x ) = KV'ilV'tT) 9))\ 2 . Then the error probability corresponding 
to the value 9 is: 

Perr(T,9) = 1 - \(^(T, 9))\ 2 

1 Quantum measure fundamental postulates can be founded in standard quantum me- 
chanics textbooks, see e.g. or 0. 
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= (^{9,T)\^{9,T)), (3) 



thanks to the fact that \i[)(T,9)) is normalized. By assumption, we 
have P err (T,9*) = 0. 

2 Robustness of the control strategy 

In the quantum control field, the expression "robustness of the con- 
trol strategy" means that the control performance is insensible with 
respect to errors in the control implementation. In [7j, a control strat- 
egy is considered robust u if significant local changes in the amplitude 
and the form of the pulse and of the chirp do not change significantly 
the final transfer probability." The pulse and the chirp, in the setting 
described in Sectional are the system inputs parameters. A quanti- 
tative definition of robustness is, however, missing. 

The contribution of the present paper is to provide such a quantita- 
tive definition reformulating the problem as a robust control-theoretic 
problem, and then to analyze the robustness properties of some sig- 
nificant strategies considered in the relevant literature. In terms of 
our model, this concept of robustness can be qualitatively formulated 
as follows: A control strategy is robust when, for values of the pa- 
rameters 9 different from the nominal ones, the final state \ip(T,9)) 
is close to the desired one \tpi). This robustness request is satisfied if 
P err (T, 9) is small in the parameter set T. 

In classical control theory, plant uncertainty is described by a set 
V of possible plants 3 . This uncertainty can be either structured 
(parametrized by a finite number of scalar parameters or a discrete 
set of plants) or unstructured (disc-like uncertainty). A controller is 
said to be robust with respect to some property if this property holds 
for every plant in V . It is quite simple to reformulate our problem as a 
particular case of structured- like classical robustness problem. First, 
we notice that our quantum "plant" is determined by the matrices 
(Hq, Hi, H m ). These matrices determine the system Hamiltonian 
(J2J), given the control strategy. Let Pq = (Hq, H%, H m ) be our 
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nominal plant. 

As in |Hj, we can transfer the uncertainty from the control pa- 
rameters to the internal Hamiltonian. In fact, by defining 8ui{9) = 
Ui(9) — Ui(9*) we can write: 

m 

H(t) = # o +^#iM0*) + M0)) 
»=i 

m m 

= H + J2 H iSui(9) + J2 H M&*) 

i=i 1=1 

m 

= (H + AH u (9)) + ^H iUi (9*). (4) 
i=i 

where AH U (9) = ^2d = iHi5ui{9). Such a cosmetic transformation 
shows that our control strategy uncertainty can be seen as a par- 
ticular case of the plant uncertainty (with control inputs Ui{9*)). The 
plant set V is here given by: 

V = {(H + AH U (9),H U ...,H m )\9 e T}. 

The property we are interested in, as e.g. in [Jj, is the error probability. 
We require this probability not to exceed a fixed threshold e G [0, 1) at 
a given T. All the ingredients of a classical robustness problem have 
now been specified. Introduce the e-robustness set 1Z t as 

K e = {9£T\P err (9,T)<e}. (5) 

We give the following definition. 

Definition 1. A control strategy {u\{9*, t), u m (9*, t)}, t £ [0, T] is 
e-robust with respect to parameters uncertainty if: 

K e = T. (6) 

Notice that only the 0-robustness case ensures us an exact steering of 
the system state for all 9 £ T. 
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3 Some applications 



In this section we analyze, in the light of the above definition, the 
robustness of some prototype examples. We will compare our results 
with qualitative observations and robustness claims in the relevant 
literature. To do so, we introduce a particular form of Q frequently 
used in NMR (Nuclear Magnetic Resonance) quantum control prob- 
lems. We consider a two level quantum system, and the associated 
a bi-dimensional Hilbert space. The time evolution is described by a 
scaled time Schrodinger equation in the form: 



ih^(s))=TH(sM(s)), (7) 



where s = t/T and 



H(s) 



-A(s) n(s) 
n(s) A( s ) 



is represented in the canonical (diabatic) base. The control functions: 

( A(s) = A <S>(s) 

{ n( s ) = n A( s ) ' 

are the inputs, with <I>(s),A(s) fixed envelops and Ao, 0,q £ ^ + am- 
plitude parameters. In this picture we have A(s) = ui(s,Ao) and 
£l(s) = U2(s,f2o). Thus 8 = (61,82) = (Aoj^o) are the parameters we 
are interested in. In the contest of particle-laser field interaction and 
the RWA (Rotating Wave Approximation pQ), these functions depend 
on the chirp (detuning) and the amplitude (time- dependent Rabi fre- 
quency) of the active pulse 2 . This model can be seen as a particular 
case of model (^Q), and it is suitable to describe control techniques 
based both on magnetic resonance and adiabatic passage. At each 
time the unitary transformation 

, „ / cos#(s) — sin 9(s) \ , s 

U(S) = • flM ' (8) 

\ smU(s) cost/(sj I 



2 To find some detailed information about the physical meaning of these parameters 
and about the resonance phenomenon see 0,01- 
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with 9(s) = | arctan (0(s)/A(s)), diagonalizes the Hamiltonian 

uHs)H(s)U(s) = ( £(S) ° / )=D(s). (9) 





Here 



e(s) = v / A2( s ) + 02 (s) . (10) 

is the energy eigenvalue. Applying U(s) as a time dependent change 
of basis and defining </>, we obtain: 

.9|0(s)) f„„, , , , d 



\<t>(s)) 
a)). (11) 



/ Te(s) i 7 ( a ) \ 

V i7(«) J 
The new basis vectors are called adiabatic states. In the adiabatic 
limit, T — ► oo 3 the 7(5) terms can be neglected, as shown by the 
adiabatic approximation theory 0. 

In the following subsections, we will investigate the robustness 
properties of different control strategies in a typical steering problem, 
the state- flipping: Transfer the system state from one basis vector 
to the other. The standard resonance technique and two adiabatic 
models will be discussed and compared. 



3.1 Magnetic Resonance 



A simple way to obtain such a transfer is by using the magnetic reso- 
nance phenomena: Under properly tailored oscillating fields, the state 
vectors rotate between the two basis states This kind of effect 

can be generated by the following fields-control functions: 

f A( S ) = 

\ n( s ) = n A( s ), 

where A(s) is the O-pulse envelope. This parametrization, and some 
easy calculations 0, lead to the following expression for the error 
probability: 



(12) 



P err (T,n ,A A ) = cos 2 



r s f 

, / A(s)ds 

J Si 



cos 2 QqTA\, 



(13) 



with A\ the 0-pulse area. 

This probability is equal to zero for: 

_ (fc + ^TT 

or, equivalently, for: 

Thus, P err (T, Slo, A\) = in a family of hyperbolas parameterized in 
k (zero- measure set in the parameters space). Consider T = [Qq — 
(3, + (3] x [A* A — a, A* A + a] , a common setting if we are working with 
nominal values and error intervals, with (3, a such that sin 2T£IqAa is 
monotone in every direction. Then, the maximum absolute value for 
the error probability in T is: 

Pmax =COS 2 TfL4 A , (14) 

where A\ = A\ + a and Q = Qq + (3. Then T is e-robust, with 

£ = Pmax- 

According to qualitative evaluation found, the magnetic resonance 
strategy doesn't seem to ensure enough insensibility towards errors in 
control implementation and can be sensitively improved by adiabatic 
passage techniques. 



3.2 Landau-Zener Model 

The Landau-Zener model is one of the simplest choice of control func- 
tion leading to an adiabatic transition. We will consider: 

(15, 

the detuning varies linearly with a zero crossing, while the Rabi fre- 
quency is maintained constant. For s = we have a minimum in 
the difference between energy levels that leads to a state inversion if 
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the evolution satisfies the condition needed for the adiabatic approx- 
imation. The error probability is estimated with the Landau-Zener 
formula: 

n 2 

Perr(T 5 Oo ! Ao)«e" 7rT 4. (16) 

This probability goes to zero in the adiabatic limit T — > oo for any 
choice of $7 7^ and . Thus, the robustness set for this strategy is 
the whole open first quadrant without its boundary (J)q = 0, Ao = 0). 
The advantages given by this adiabatic technique are evident, as long 
(j!6|) estimates correctly the error probability. Even if the T — » 00 
condition is not realizable, we can take a T large enough to maintain 
P err arbitrary small for (almost-)every parametrization of the control 
strategy. We will call this behavior intrinsically robust. 



3.3 Allen-Eberly Model 

We now analyze an adiabatic control strategy more complex than the 
previous one. The Allen-Eberly parametrization allows to obtain 
an exact expression for the error probability and, in the Qq = Aq case, 
forces the state time evolution along the energy level lines, maintaining 
e(0(s), A(s)) = c, c constant for every s [Zj. This kind of choice leads 
to good results in terms of error probability even quite far from the 
ideal T —* 00 condition, as we are going to show. In terms of control 
functions, we consider: 

A(s) = A ^ 1 - sech 2 (s) = A tanh(s) ^ 
f2(s) = f2osech(s). 

Then, the exact expression for the error probability is: 

P err {T, O , Ao) = cosh 2 ^vrT^A 2 - Q^j sech 2 (vrA T) , (18) 

for every regime, adiabatic or not. We can notice that, for large T 
and for Ao > Oo> the error probability can be bounded by: 

P err (T,n ,A ) < 4 e -2TT(Ao-VAg-ng)_ (19) 
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Thus, for every Ao and CIq, Aq > Qq, the error probability decreases 
exponentially to zero in the adiabatic limit. The best choice for the 
parameters values is to take the largest Ao = fio- I n the case £Iq > Ao, 
the error probability becomes: 



P err (T, Sl , A ) = cos 2 ( irTy fig - A 2 J sech 2 (vrA r) . (20) 
This expression tends to zero with dumped oscillations, due to the 



term cos 2 yirTy/ClQ — AqJ. Again, larger Ao make P err converge 
faster. Thus, for each fixed e, we can compute a T e such that the 
error probability P err (T, CIq, Aq) < e for every T > T e . Indeed, it is 
easy to see that 

rp r ^ lne 

l f = max J — , — \ . 

2vr(A -7A^2)' 2vrA ; 

This control strategy is therefore intrinsically robust for T sufficiently 
large. According to the Landau-Zener case, every choice of S7o 7^ and 
Ao 7^ drives the system to the target state. The level line condition 
(Aq = fio) and large Qq give faster convergence to the desired state. 



4 Discussion 

Comparing the results, the advantages given by the adiabatic strate- 
gies are evident. They can be effectively used, however, when the 
transfer time is not critical: Their intrinsic robustness is exhibited 
only with a large time. 

The examples analyzed are also treated in to illustrate that con- 
trol strategies based on the level lines are optimal for adiabatic popu- 
lation transfer (the level line strategies minimize the error probability 
and corresponds to the minimum pulse area). In [Jj, robustness of 
the control is also taken in account: Contour plots of error probability 
with respect to parameters variations are obtained thanks to numerical 
simulations for the system evolution. Different strategies are qualita- 
tively compared. It is shown that the simple resonance case generates 
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larger error probability than the adiabatic optimal techniques, once a 
parameter variation is fixed. 

Here we have given a formal definition of the robustness property, 
reformulating the problem in the control theoretical setting. We have 
obtained quantitative results consistent to the qualitative ones just 
described, and we have provided an analysis tool useful to evaluate and 
compare robustness behavior of different strategies. From a control 
theoretic viewpoint, we have analyzed a specific robustness problem 
for open-loop control of a bilinear system. 
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